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Gyrokinetic validation studies are crucial in developing confidence in the model incorporated in numerical simula¬ 
tions and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade 
discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measure¬ 
ments. The approach taken is the following. Eirst, linear simulations are performed in order to determine the turbulence 
regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic 
ion temperature gradient within the expected experimental error bars. Einally, the dependence of various quantities 
with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations 
can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive 
to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better 
observables for comparisons between gyrokinetic simulations and experimental measurements. 


I. INTRODUCTION 

It has been known for several decades that energy and par¬ 
ticle confinement in tokamak plasmas are mainly degraded 
by turbulence driven by steep temperature and density gradi¬ 
ents. Eor this reason, the characterization and understanding 
of these turbulent processes is a very important task in order 
to improve the performance of present experiments as well as 
future fusion reactors. 

Due to the strong background magnetic held and the low 
collisionality in tokamak plasmas, gyrokinetic theory has been 
established as the most appropriate theoretical framework for 
the study of turbulent transport in the plasma core. In order 
to improve confidence in the numerical results obtained with 
gyrokinetics and to establish a solid understanding of turbu¬ 
lent transport across the whole range of plasma parameters, 
it is very important to perform direct comparisons between 
simulations and experimental measurements. In this respect, 
with the recent development and improvements in fluctuations 
diagnostic, it is now possible to measure turbulence features 
with high precision, allowing for quantitative comparisons be¬ 
tween experimental data and results of nonlinear gyrokinetic 
simulations. These validation studies are crucial in developing 
confidence in the models and improving the predictive capa¬ 
bilities of the numerical simulations. 

In a recent paper—, we have compared density fluctuation 
levels measured with a new Doppler reflectometer installed 
in ASDEX Upgrade and simulation results obtained with the 
gyrokinetic GENE code. We extend the previous work by ad¬ 
ditionally presenting simulation results of density wavenum¬ 
ber spectra, electron temperature fluctuation levels as well as 
cross-phases between different quantities. One of the rea¬ 
sons for analyzing electron temperature fluctuation levels and 
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cross-phases is that a new Correlation Electron Cyclotron 
Emission (CECE) system is expected to be installed and to be 
in operation in 2015 in ASDEX Upgrade. Therefore, the gy¬ 
rokinetic results presented in this paper can provide guidance 
for the on-going development of the diagnostic. Erom a more 
fundamental point of view, we will also investigate the varia¬ 
tion of these quantities with respect to various physical input 
parameters. This information can be used to characterize core 
turbulence features in ASDEX Upgrade plasmas. 

The paper is organized as follows. In Section [III an 
overview of the chosen plasma discharge analyzed is given. 
A description of the gyrokinetic simulation method used is 
described in Section |III1 Micro-instability studies from lin¬ 
ear gyrokinetic simulations are outlined in Section HV] The 
main results of the paper are shown in Section|V] Core turbu¬ 
lence features such as heat fluxes, density fluctuation ampli¬ 
tudes and spectra, temperature fluctuation amplitudes as well 
as cross-phases between these quantities will be presented in 
detail, followed by a discussion in Section |VT] Einally, con¬ 
clusions and future work will be discussed in Section lVlIl 


II. OVERVIEW OF THE PLASMA DISCHARGE 

The ASDEX Upgrade discharge ^28245 analyzed in this 
paper was operated in the high-confinement regime (H-mode). 
It was planned to study the turbulence characteristics in both 
the ion temperature gradient (ITG) and trapped electron mode 
(TEM) regimes, through a transition from one regime to an¬ 
other. This transition can be achieved by modifying locally the 
electron temperature gradient, which affects the TEM instabil¬ 
ity. This is obtained by changing the electron cyclotron reso¬ 
nance heating power (ECRH) by steps during the discharge. 

An overview of several relevant time traces is given in 
Fig.lU The ECRH power Pecrh is deposited at ppoi = 0.5 
from 2.0 — 4.0 seconds. Here, ppoi is the normalized poloidal 
flux radius. At this heating location, Pecrh is varied sub- 
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TABLE I. Physical parameters for the six simulated ASDEX Up¬ 
grade cases. i?axis is the major radius at the magnetic axis and 
Ps = Cs/0,i is a reference gyroradius defined with the ion sound 
speed Cg = \/Tglmi and the ion gyrofrequency f2i. 
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FIG. 1. (color online) Time traces of ASDEX Upgrade discharge 
^28245. (a) ECRH heating power, (b) electron temperature and 
(c) logarithmic electron temperature gradient at different radial posi¬ 
tions. The data is analyzed in the time intervals shaded in grey. 


sequently between 0.5,0.0,1.2 and 1.8 MW (a). The influ¬ 
ence of the stepped heating power can be clearly observed 
in the electron temperature (b). Finally, the strongest in¬ 
crease of the logarithmic temperature gradient uit^ is observed 
at Ppoi ~ 0.6 when an ECRH power of 1.8 MW is applied 
(c). For this reason, this case and the one without ECRH will 
be analyzed in detail, as examples of the two extreme cases. 
They correspond to the time windows shaded in grey in Fig.[T] 
Within these time windows, we have simulated three different 
radial locations: ppoi = {0.6,0.7,0.8}, making a total of six 
different scenarios to be studied with gyrokinetic simulations. 


III. OVERVIEW OF THE GYROKINETIC SIMULATION 
METHOD 

The turbulence data obtained in this paper are pro¬ 
duced with the gyrokinetic code GENE, which solves self- 
consistently the gyrokinetic-Maxwell system of equations on 
a fixed grid in five dimensional phase space (plus time): two 
velocity coordinates (u||, /r) and three field-aligned coordi¬ 
nates (x, y, z). Here, z is the coordinate along the magnetic 
held line, while the radial coordinate x and the binormal coor¬ 
dinate y are orthogonal to the equilibrium magnetic held. The 
velocity coordinates are, respectively, the velocity parallel to 
the magnetic held and the magnetic moment. GENE has the 
possibility to simulate either a flux-tube (local simulations) or 
a full torus (global simulations). In the former option, it is 
assumed that the relevant turbulent structures are small with 
respect to the radial variation of the background profiles and 
gradients. This allows to use periodic boundary conditions 
and thus the coordinates perpendicular to the magnetic held 
are Fourier transformed (a:, y) —)• {kx,ky). For this work, 
only the local version of the code has been employed. 


The physical parameters for each of these six cases are 
given in Table U There, the reference length is defined as 
Tret = \/'E'tor.sep/7ri?ref, where fitter,sep is the toroidal flux 
at the separatrix and Sref is the magnetic field on axis. Typ¬ 
ically, this reference length is comparable but not identical to 
the tokamak minor radius. The logarithmic gradients are de¬ 
fined as in Ref. H: ojx = ^ ^ {Ti,Te,ne}, 

and ptor the normalized toroidal flux radius. The magnetic 
shear is given as s = where q is the safety factor, 

and the electron beta is defined as /3e = 2/iotteTe/T^ef' Here, 
Ti and Tg are the ion and electron temperature, respectively, 
and Tie is the electron density. 


The GENE code is physically quite comprehensive and in¬ 
cludes many features (see Ref. for more details). For the 
ASDEX Upgrade scenario studied in this paper, the following 
features of the GENE code were used: two particle species 
(deuterons and electrons), electromagnetic effects by solving 
the parallel component of Ampere’s law, external ExB shear, 
parallel flow shear and a linearized Landau-Boltzmann colli¬ 
sion operator with energy and momentum conserving termsi. 
Unless stated otherwise, the magnetic equilibrium geometry 
is taken from the TRACER-EFIT interface^. Additionally, 
GyroLES techniques have been used to reduce the accumula¬ 
tion of energy at the smallest scales (see Refs. |@-[^). Further 
simulation details, such as resolution grid, box sizes, etc., are 
given in the following sections. 
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IV. MICRO-INSTABILITY STUDIES: LINEAR 
GYROKINETIC SIMULATIONS 

In order to calculate turbulent transport fluxes, density 
and temperature fluctuation amplitudes, etc., nonlinear gy- 
rokinetic simulations are necessary. Nevertheless, linear gy- 
rokinetic simulations can provide useful insights. For in¬ 
stance, they may allow us to identify the underlying micro¬ 
instabilities which drive the turbulence present in the exper¬ 
iments. They can also be used for convergence studies and, 
since they are usually computationally cheap, they can also be 
used to do scans in different physical parameters. 

1. Nominal parameter set 

In linear simulations we calculate the growth rate and fre¬ 
quency of the most unstable mode present in the system for a 
given binormal wave vector ky and kx = 0. In this paper we 
choose to present ky in cm~^ instead of the more common 
kyPs units. This has been done in order to compare to experi¬ 
mental results. For reference, ps values are given in TableU 

In Fig |2] we display the growth rates ( 7 ) and frequencies 
(w) for the cases at ppoi = 0 . 6 , because similar conclusions 
are obtained for the other cases. In the figure, the negative fre¬ 
quencies are represented by dashed lines. The grid resolution 
was {x, z, v\\,p} = {31, 32,48,16}. Convergence tests were 
performed at higher resolutions and confirm the validity of 
the results. Several observations can be made. First, for low 
wavenumbers, ITG is the dominant instability. This is indi¬ 
cated by a positive frequency, which with the present normal¬ 
ization represents a frequency in the ion diamagnetic direc¬ 
tion. Second, for all the cases analyzed, the TEM mode is sta¬ 
ble (studied with an eigenvalue solver). In fact, as was already 
shown in Ref. [I], only with the combination of a much higher 
electron temperature gradient and a lower ion temperature gra¬ 
dient than the ones measured experimentally, do TEM modes 
become unstable^. Third, for higher wavenumbers, ETG is the 
dominant instability (indicated by a negative frequency). Ei- 
nally, there is practically no effect of the ECRH on the ITG 
growth rates. The main effect of ECRH is to increase the 
growth rates of ETG modes, possibly leading to a subsequent 
increase of the electron heat flux for these cases. However, in 
both cases, we expect their contribution to be small with re¬ 
spect to the ITG contribution. Moreover, ETG modes are not 
expected to influence density and temperature fluctuations at 
low wavenumbers, which are dominated by ITG. Since these 
are the scales measured by the diagnostics we are considering 
here, in this work we will limit ourselves to wavenumbers up 
to ky = 10cm“^, thus excluding ETG modes, but allowing 
for a significant reduction in computational resources. 

2. Sensitivity studies with respect to the main physical 
parameters 

Physical parameters such as temperature, density, magnetic 
equilibrium profiles, etc., are measured with experimental un- 



FIG. 2. (color online) Linear growth rates (left) and frequencies 
(right) for the cases at ppoi=o .7 versus the binormal wavenumber 
ky. The dashed lines in the figure indicate negative frequency values. 
Positive frequencies refer to modes drifting in the ion diamagnetic di¬ 
rection and negative in the electron diamagnetic direction. For these 
cases, they correspond to ITG and ETG modes, respectively. 


certainties. Since these values are used as input in the gyroki- 
netic codes, the uncertainty in these quantities could affect the 
simulation results. Therefore, sensitivity studies are carried 
out, with the aim of studying the effect of the different uncer¬ 
tainties on the simulation results. Ideally, these studies should 
be done in nonlinear gyrokinetic simulations. However, due to 
the expensive computational effort associated, this is in prac¬ 
tice unfeasible. For this reason, this sensitivity study is done 
within linear gyrokinetic simulations. 

We have studied the sensitivity of the linear growth rate 
with respect to a variation of ±20 % in the nominal value 
for different physical parameters, such as: logarithmic ion 
temperature gradient (wtJ, logarithmic electron temperature 
gradient (uiTe), logarithmic electron density gradient (w„,,), 
electron to ion temperature ratio {T^/Ti), collisionality (r'coi), 
safety factor (q) and magnetic shear (s). The main results are 
summarized in Table [III For simplicity, we show only cases 
at Ppoi = 0.7, although similar conclusions were obtained to 
the other radial positions. The sensitivity studies for q and 
for s have been done using a Miller-type magnetic equilib- 
riumiS. As is shown in Table HIl the peak of the growth rates 
are practically insensitive with respect to changes of ± 20 % in 
ujn, ^Te, i^coi, 9 and s. Changes in the peak of the growth rate 
up to 20 % are found for T^jTi variations. The most critical 
parameter is whose ± 20 % variation modifies the growth 
rate by up to 40% (see Fig. (3). Based on these results, one 
could expect that the uncertainties in iot^ will have the largest 
influence in nonlinear gyrokinetic simulations. In the follow¬ 
ing, we will mainly focus on the influence of this parameter 
on nonlinear simulations. The influence of the E x B shear 
is in general also expected to have a relevant impact on the 
transport in nonlinear simulations. However, the low value of 
the E X B shear for this particular discharge is such that its 
influence can safely be neglected. 
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No ECRH 

A 7 

1.8 MW ECRH 

A 7 

OJTi X 0.8 

-38 % 

OJTi X 0.8 

-30% 

UJTi X 1.2 

+25 % 

UJTi X 1.2 

-1-26% 

LJTe X 0.8 

-4% 

iOTe X 0.8 

-5% 

iOTf, X 1.2 

-1-5 % 

X 1.2 

-e5% 

Wrig X 0.8 

-1-0% 

Wrig X 0.8 

-5% 

OJn^ X 1.2 

-1-0% 

^Tie X 1.2 

-1-6% 

T^/Ti X 0.8 -1-14% 

Te/Ti X 0.8 

-Ell % 

Te/Ti X 1.2 

-18 % 

Te/Ti X 1.2 

-9% 

g X 0.8 

-14% 

g X 0.8 

-9% 

g X 1.2 

-1-5 % 

g X 1.2 

-e6% 

s X 0.8 

-1-0% 

s X 0.8 

-5% 

a X 1.2 

-4% 

s X 1.2 

-e5 % 

i^coi X 0.8 

-1-5 % 

^^coi X 0.8 

-e6% 

i^coi X 1.2 

-4% 

^coi X 1.2 

-5 % 


TABLE II. Percentage difference in maximum growth rate with re¬ 
spect to the nominal values for cases at Ppoi = 0.7 for various pa¬ 
rameters (curi, uJTe, uin, TijTe, g, s, and r'coi)- 


no ECRH 1.8 MW ECRH 



FIG. 3. (color online) Linear growth rates versus the binormal 
wavenumber ky for the cases at ppoi = 0.7 with respect to the varia¬ 
tion of the logarithmic ion temperature gradient. 


V. CORE TURBULENCE FEATURES: NONLINEAR 
GYROKINETIC SIMULATIONS 


In order to predict and compare with experimental results, 
nonlinear simulations are required. For the selected discharge, 
the grid resolution needed is {256 x 128 x 32 x 48 x 16} 
points in {x, y, z, un, p} coordinates. A convergence test of 
the results with that resolution has been performed by com¬ 
parison with nonlinear simulations with a double resolution in 
the perpendicular directions for a few cases. Moreover, per¬ 
pendicular box sizes have been chosen in such a way that sev¬ 
eral correlation lengths fit in the box and convergence checks 
have also been done on this respect to ensure the validity of 
this choice. Results of the simulations are time-averaged over 
a range well exceeding the correlation time of the underlying 
turbulence. 


A. Turbulence ion and electron heat fluxes 

In this section, we compare the experimental ion and elec¬ 
tron heat fluxes obtained through power balance analysis with 
the ASTRA code and the results from nonlinear GENE simu¬ 
lations. As was done previously, we have grouped in the same 
plot the cases taken at the same radial position. The results 
are shown in Eig. ID 
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FIG. 4. (color online) Comparison of experimental heat fluxes 
(dashed-lines) with those obtained from gyrokinetic nonlinear GENE 
simulations (markers) using the nominal parameters and the variation 
with respect to logarithmic ion temperature gradient for all cases. Ion 
heat fluxes (left) and electron heat fluxes (right) are depicted. The 
rows represent the radial positions ppoi = 0.6, 0.7, 0.8. The experi¬ 
mental values are obtained through power balance analysis with the 
ASTRA code and the shaded regions are used to indicate the uncer¬ 
tainty of the ASTRA values. In this case, a 20% error is assumed for 
all cases. In blue are the discharge parameters without ECRH heat¬ 
ing and in red with ECRH heating. The nominal parameter in each 
case is colored differently to distinguish it from the rest. 

Ion (electron) heat fluxes are shown in the left (right) 
columns in Eig. 01 The rows represent radial positions ppoi = 
0.6, 0.7,0.8. The dashed lines indicate the ASTRA results and 
their shaded regions are used to indicate the uncertainty of the 
ASTRA values, where a 20 % uncertainty is assumed^. Based 
on the linear sensitivity studies, for each case, several simula¬ 
tions were performed varying the logarithmic ion temperature 
gradient in steps of ±10% with respect to the nominal values, 
up to a maximum of ±30% variation. GENE simulation re¬ 
sults are represented by the markers in the figure. The cases 
without ECRH are colored in blue and with ECRH in red and 
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for each set, the simulations with the nominal parameters are 
colored differently (in green for the cases without ECRH and 
in yellow for the cases with ECRH). The statistical error bar is 
an estimation of the standard deviation of the set of means of 
consecutive temporal sub-domains of the saturated state. Sev¬ 
eral conclusions can be obtained from Eig.|4] Eor the nominal 
parameters, the cases with ECRH produce more ion and elec¬ 
tron heat flux than the ones without ECRH for all positions. 
Moreover, at ppoi = 0.6, the ion heat fluxes match the experi¬ 
mental values without having to vary the gradient with respect 
to the nominal value. At this position, only the electron tem¬ 
perature must be increased by 10%. However, for the other 
positions, the values of the heat fluxes obtained with the nom¬ 
inal parameters clearly overestimate the heat fluxes obtained 
with ASTRA by a factor of 2 — 3. We need to decrease wy. by 
20 % in order to match the experimental results for the cases 
without ECRH. Whilst, for the cases with ECRH heating, it 
has to be decreased by a maximum of 30%. In Eig. |5] only 
the flux-matched results are compared with the experimental 
heat fluxes. We can conclude that agreement of the transport 
levels within the errors bars can be achieved, since even small 
uncertainties in the temperature profile itself may translate to 
relatively large ones (up to 20 — 30%) in the logarithmic gra- 
dientsiii^. 



/^pol /^pol 

FIG. 5. (color online) Comparison of experimental heat fluxes 
(dashed-lines) with those obtained from gyrokinetic nonlinear GENE 
simulations (markers) using the flux-matched simulations. 


1. Turbulence electron density amplitudes 

Since turbulent fluxes are caused by plasma fluctuations 
on microscopic scales, it is necessary to validate gyrokinetic 
codes on a microscopic level. In this regard, a new Doppler 
reflectometer has been recently installed in ASDEX Upgrade 
(see Ref[I]for more details on the diagnostic), which is able to 
measure electron density fluctuation amplitudes (fl-e). In order 
to compare experimental and numerical results, a synthetic di¬ 
agnostic must be implemented in GENE to reproduce the mea¬ 
surement process of the reflectometer. Two kinds of synthetic 
diagnostics can be employed. A first approach consists in sim¬ 
ply Altering the data in order to taking into account only the 
location and the wavenumbers that the diagnostic measures. 
This translates to take into account only fluctuations at the out¬ 
board mid-plane (z = 0), averaged over a finite radial length 
and then to select the range of perpendicular wavenumber that 
are measured in each case < ky < 


GENE). A more sophisticated method uses a full-wave code 
to simulate also the incidence and reflection of the wave into 
the gyrokinetic turbulent dataii. This work is in progress and 
only preliminary results are available with this synthetic diag¬ 
nostic (see Ref. [fsl) . Eor this reason, we have only used the 
Altering method for this work. 

Comparison of the experimental and simulated fie are 
shown in Eig. |6l The simulations that match the experi¬ 
mental ion heat fluxes are shown with a different marker to 
distinguish them from the rest. The fluctuation data is ana¬ 
lyzed considering only perpendicular wavenumbers between 
A < ky <8 [cm“^]. The experimental values are scaled by a 
common factor since the measurements are in arbitrary units. 
Eor this reason, only the shape of radial turbulence level pro¬ 
files and the effect of ECRH can be used for comparison. We 
decided to scale the experimental values to try to match the 
case without ECRH. As is shown in the left plot of Eig. |6] 
we obtain a remarkable agreement between experimental and 
simulations results in the radial trend. Eor the case of 1.8 MW 
ECRH, there is also a good agreement in the turbulence level 
profile. However, with the scale used, the fluctuation levels 
are clearly underestimated with respect to the ones measured 
experimentally. In particular, the flux-matched results present 
the biggest discrepancy with respect to the experimental mea¬ 
surements. Einally, from this figure we can also observe how 
sensitive the density fluctuations are with respect to variations 
in the ion logarithmic temperature gradient. Eor instance, a 
30% reduction in the logarithmic gradient can reduce density 
fluctuation levels by more than a factor of 2. 
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FIG. 6. (color online) Electron density fluctuation amplitudes at dif¬ 
ferent radial positions. Blue-plus markers represent the data without 
ECRH heating, red-cross markers with ECRH and experimental re¬ 
sults are in full circles. The flux-matched (GENE f.m.) simulations 
are marked differently. Eluctuation data is analyzed at the outboard 
mid-plane, averaged over a finite radial length and with perpendicu¬ 
lar wavenumbers between 4 < ky <8 [cm“^]. 


2. Turbulence electron density spectra 

The knowledge of the power-law spectra of a physical 
quantity is important for the understanding of the underlying 
physics and useful for providing constraints for simple phys¬ 
ical models. Based on Kolmogorov-type argumentsi^, turbu¬ 
lence is generally associated with universal power-law spec- 
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FIG. 7. (color online) Electron density fluctuations spectra at ppoi = 
0.8 including the variation of the logarithmic ion temperature gradi¬ 
ent. The solid lines represent the wavenumber range where the fit to 
a power law was done. Spectral indices for the flattest and steepest 
spectra are also indicated in the figure for each case. 



ijjrp Ccl'J' 


FIG. 8. (color online) Calculated spectral indices versus the logarith¬ 
mic ion temperature gradient for different radial positions and ECRH 
scenarios. A decrease of the spectral indices with respect to the ion 
temperature gradient is observed for most of the cases. 



FIG. 9. (color online) Calculated spectral indices versus the radial 
position. Blue-plus markers represent the data without ECRH heat¬ 
ing and red-cross markers with ECRH. The flux-match simulations 
are marked differently. A decrease of the spectral indices with re¬ 
spect to the radial position is observed. 

tra. However, as shown in Refs and this is gener¬ 

ally not the case in plasma turbulence, and different power- 
law spectra indices can be found depending on the type of 
mechanism which drives or dissipates energy in the system. 
Furthermore, the knowledge of wavenumber spectra could be 
important for a clear identification of the turbulent regimes 
driven by different microinstabilities and can be used to fur¬ 
ther validation of the gyrokinetic model. 


The results of the electron density fluctuation spectra at 
Ppoi = 0.8 for different logarithmic ion temperature gradients 
are shown in Fig.]?] The solid lines represent the wavenumber 
range where a fit to a power law is shown: (|neP) = ak^, 
where b is the spectral index, and () represents an average over 
all the coordinates except k±. Spectral indices for the flattest 
and steepest spectra are indicated in the figure. The calculated 
spectral indices are also shown in Figs. 0 and |9]with respect 
to the ion temperature gradient and to the radial position, re¬ 
spectively. These figures present a clear qualitative behavior: 
the spectral index decreases with the increase of turbulence 
drive. Additionally, the spectral indices also decrease when 
going from the inner to the outer core position. Moreover, 
the magnitude of the exponents cover a wide range of values, 
approximately from 4 to 9. Although density fluctuation spec¬ 
tra were not measured with the Doppler reflectometer for this 
discharge, a similar trend has been also reported in Ref, fl^ for 
various ASDEX Upgrade discharges. 

These results indicate that the turbulence driven by ITG 
modes exhibits non-universal power laws, whose spectral in¬ 
dices could depend on several physical parameters. Future 
work in this respect will be to study also if TEM modes ex¬ 
hibit similar properties. If this was the case, then it would 
become very difficult to distinguish a type of instability by 
measuring only its characteristic spectral index. 


B. Turbulence electron temperature amplitudes 

At the time the discharge was performed, no temperature 
fluctuation measurements were available. However, a Corre¬ 
lation Electron Cyclotron Emission (CECE) diagnostic is cur¬ 
rently installed on ASDEX Upgrade, and electron temperature 
fluctuation profiles will be available in the future campaigns. 

The CECE diagnostic measures perpendicular electron 
temperature fluctuations {T±^e) in the long wavelength range 
(relevant for ITG and TEM modes) and is not sensitive short 
wavelengths (ETG modes). This diagnostic presents an inher¬ 
ent limitation in the lowest fluctuation level that can be de¬ 
tected. This noise level depends on the physical parameters of 
the discharge, but typical values are between 0.2 — 0.3 %2a. 

A detailed description of CECE modeling in DIII-D is 
given in Ref. This synthetic diagnostic has already been 
implemented in GENE for DIII-D discharges'^ and a simi¬ 
lar synthetic diagnostic will be implemented for ASDEX Up¬ 
grade discharges. However, this diagnostic could be not used 
in this work since it requires the knowledge of the CECE con¬ 
figuration during the discharge. Eor this reason, we have con¬ 
sidered a simpler synthetic diagnostic, which consists in fil¬ 
tering the gyrokinetic data to the positions and wavenumbers 
that are expected to be measured in ASDEX Upgrade. Con¬ 
sequently, the gyrokinetic data analysis results are restricted 
to the outboard mid-plane position {z = 0), averaged over the 
finite radial length and summing all perpendicular wavenum¬ 
bers (since short wavelengths are not simulated). In order to 
better model the actual diagnostic, one should also apply a 
filter in the frequency space. However, this filter will also de¬ 
pend on the specific range of frequencies measured. Since we 
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do not have access to this information, we have considered all 
the frequencies in the analysis. Therefore, the following re¬ 
sults should be only used as an approximated indication of the 
fluctuation amplitudes that could be detected with CECE for 
this discharge. Nevertheless, we do not expect radial trends 
to change with respect to a more sophisticated synthetic di¬ 
agnostic approach and only the amplitudes are likely to be 
rescaled^^. 

The main results are shown in Eig. [TOl The perpendicular 
electron temperature fluctuation amplitudes go from a mini¬ 
mum of 0.1 % at the inner position to a maximum of 1.0 % at 
the outer core position. Therefore, assuming a noise-level in 
the order of 0.2 — 0.3 % for the CECE diagnostic, this could 
imply that only the fluctuations at the outer core positions 
(starting from ppoi > 0.6) could be detected. As for the case 
of the density amplitudes, we also observe a large variation of 
the fluctuation amplitudes with respect to the changes in the 
logarithmic ion temperature gradient. 



FIG. 10. (color online) Percentages of electron perpendicular tem¬ 
perature fluctuation amplitudes at different radial positions. Blue- 
plus markers represent the data without ECRH heating and red-cross 
markers with ECRH. The flux-match simulations are marked differ¬ 
ently. The area shaded in grey indicates a typical noise level of the 
CECE diagnostic. 


C. Turbulence cross-phases 

Doppler reflectometers can be coupled to the CECE diag¬ 
nostics to calculate cross-phases between electron density and 
temperature fluctuations^. This measurement is important for 
gyrokinetic validation studies since it represents a relationship 
between different fluctuating quantities (density and tempera¬ 
ture in this case). In addition, this cross-phase could be also 
related to the cross-phase that determines the turbulent heat 
fluxes (electrostatic potential and temperature fluctuations). 
Eor this reason, in addition to the cross-phase that can be 
measured experimentally, we will also show the cross-phase 
between electrostatic potential and electron density fluctua¬ 
tions, so we can relate the cross-phases measured experimen¬ 
tally to the turbulent heat fluxes. The cross-phases are here 
defined as (5-^ = ta,n~^{^{h/Tj_^e)/^{n/Tj_^e)), and, 
= tan-i(Ss(n/^)/sR(n/(^)). 

Eig. [TT] shows the variation of the cross-phase versus the 
ion temperature gradient integrated over binormal wavenum¬ 
bers in the range of 1 cm“^ < ky < 10cm“^. Eor most of 



UJ'j’ LO'j’ 


FIG. 11. (color online) Calculated cross-phases for the different 
cases versus the logarithmic ion temperature gradient. Cross-phases 
between density and electrostatic potential fluctuations (a) and (b). 
Cross-phases between electron density and temperature fluctuations 
(c) and (d). The cross-phases seem to be rather insensitive with re¬ 
spect to changes in the logarithmic ion temperature gradients. 


the cases, the phases remain rather invariant with respect to 
the variation of this parameter. Therefore, this result seems to 
indicate that the cross-phase is a better observable to identify 
the type of instability which drives the turbulence in experi¬ 
ments, since TEM instability is expected to exhibit different 
cross-phases^. 

In Eig. [12] the cross-phase are displayed versus the radial 
positions. Regarding the cross-phase between density and 
electrostatic potential fluctuations (a) and (b), we see that they 
are practically in phase, i.e. close to 0. In addition, an in¬ 
crease of the phase with the radial position is also observed, 
going from practically 0 degrees at ppoi = 0.6 to approxi¬ 
mate by 25 degrees at the outer core position. On the contrary, 
for the cross-phase between electron density and temperature 
fluctuations (c) and (d), we see a decrease with the radial po¬ 
sition, going from around 150 degrees in the inner position 
to a value of 90 degrees, which result in an increase of the 
electron heat flux. Similar values of this cross-phase have 
also been measured in DIII-D and calculated with GYRO in 
Ref. im Eurthermore, these values have also been found by 
GENE for these discharges, see Ref. [o] These observations 
could be explained in the following way. In the inner position, 
the population of the trapped particles that contributes to the 
ITG instability is small, so the electrons behave almost adi- 
abatically. Because of this, density and potential fluctuations 
are in phase and density and perpendicular temperature fluctu¬ 
ations are almost out of phase (i.e., close to 90 degrees), thus 
producing negligible electron heat flux. With increasing the 
radial position, the population of trapped particles increase 
and a deviation of the electron adiabaticity is observed. Eor 
this reasons, both cross-phases approach to 90 degrees, with 
the subsequent increase of electron heat flux. 

Einally, in Eig. [T^ the colored contours display the cross¬ 
phases obtained from the nonlinear simulations, while the red 
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squares are used to display the cross-phases of the linear sim¬ 
ulations for the case without ECRH at ppoi = 0.6. The agree¬ 
ment between linear and nonlinear cross-phases is remarkably 
good. This is also observed for the other case a^'^^’^'* . This re¬ 
sult implies that linear simulations could be enough to com¬ 
pare with experimental results. 
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FIG. 14. (color online) Impact of the variation of the logarithmic 
ion temperature gradient around the nominal value for the case at 
Ppoi=o .7 with 1.8 MW ECRH on a) ion and electron heat fluxes b) 
electron density and temperature fluctuation amplitudes, c) electron 
density spectral index and d) cross phases. Figures a), b) and c) are 
normalized with respect to the value at the nominal ion temperature 
gradient. 


FIG. 12. (color online) Calculated cross-phases versus the radial 
position. Cross-phases between density and electrostatic potential 
fluctuations (a) and (b). Cross-phases between electron density and 
temperature fluctuations (c) and (d). The flux-match simulations are 
marked differently. 


no ECRH no ECRH 




FIG. 13. (color online) Case without ECRH at ppoi = 0.6: Com¬ 
parison of linear (markers) and nonlinear (contour) cross-phases as 
function of the binormal wavenumber where amplitudes increase 
from white to black. Left. Cross-phases between electron density 
and electrostatic potential fluctuations. Right. Cross-phases between 
density and temperature fluctuations. 


VI. DISCUSSION 

Focusing on the case at ppoi = 0.7 with 1.8 MW ECRH, 
the key results obtained in this paper can be illustrated in 
Fig- [HI Here, we show the impact of the variation of the log¬ 
arithmic ion temperature gradient around the nominal value 
on various observables. Ion and electron heat fluxes (a), elec¬ 
tron and temperature fluctuation amplitudes (b) and electron 
density spectral indices (c) are all very sensitive with respect 


to small changes in the ion logarithmic temperature gradient. 
For instance, by decreasing ujXi by 30%, the amplitudes can 
be reduced by a factor of 2 for the density and temperature 
fluctuations and by a factor of 4 for the heat fluxes. This result 
implies that the comparison of gyrokinetic simulation and ex¬ 
perimental measurements for these observables are very sensi¬ 
tive to uncertainties in the experimental input profiles. On the 
contrary, cross-phases between density and temperature fluc¬ 
tuations and between density and electrostatic potential (d) are 
rather insensitive with respect to oJXi ■ This, together with the 
fact that linear and nonlinear cross-phases agreed also remark¬ 
ably well, indicates that cross-phases could be a good observ¬ 
able to compare (fast) linear gyrokinetic simulations with ex¬ 
perimental measurements. 


VII. CONCLUSIONS AND FUTURE WORK 

We have analyzed, by means of gyrokinetic simulations 
with GENE, core turbulence features of an H-mode discharge 
in ASDEX Upgrade. The main results of this paper can 
be summarized as follows. Flux-matched simulations were 
achieved by varying the nominal ion temperature gradient by 
a factor of 20 — 30%, which is within the uncertainty range 
of the experimental profiles. In addition, density fluctuation 
levels show an agreement in the shape of the radial turbu¬ 
lence level profiles, although the effect of the ECRH on the 
fluctuation levels was not reproduced. Non-universal power- 
law spectra were found for turbulence driven by ITG modes. 
In particular, ITG instability exhibits spectral indices for the 
density fluctuation spectra which cover a broad range of val¬ 
ues. These values depend on the radial position and also on 
the specific ion temperature gradient. Gyrokinetic simulations 






































9 


predict a decrease of the exponents with respect to both the 
increase of the ion temperature gradient and the increase of 
the radial position. These results could help validate future 
analytical theories and are useful for comparisons with other 
gyrokinetic codes and future measurements. Regarding the 
electron temperature fluctuations, we observe for the inner po¬ 
sition (ppoi = 0.6) fluctuation amplitudes which are close to 
the sensitivity of the CECE diagnostic. Therefore, it is possi¬ 
ble that only measurements at positions larger than ppoi = 0.6 
could be detected with this diagnostic in such discharges. We 
hope that these results can provide guidance for the develop¬ 
ment of the CECE diagnostic that is currently being installed 
in ASDEX Upgrade. Einally, by analyzing cross-phases be¬ 
tween density and temperature fluctuations and between den¬ 
sity and electrostatic potential, we observed that linear and 
nonlinear cross-phases agree remarkably well, and that they 
are rather insensitive with respect to the variation of the ion 
temperature gradient, indicating that cross-phases could be a 
good observable with experimental measurements for com¬ 
parisons. 

Eor future work, GENE and ASDEX Upgrade comparisons 
will continue with the study of similar H-modes plasmas but 
with higher ECRH power (up to 3.6 MW). These discharges 
are expected to have peaked electron temperature profiles and 
allow TEM modes to be dominant. This will allow us to 
study fundamental differences between ITG and TEM modes 
from a microscopic level. Eurthermore, additionally dedicated 
discharges have already been conducted in which detailed 
wavenumber spectra have been measured with the Doppler 
reflectometer. These comparisons, along with the inclusion 
of future CECE measurements, will help in further validating 
gyrokinetic codes and the development of synthetic diagnos¬ 
tics. 
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